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This contract was to astablish parformanca raquiramants for 
achiaving a gaoid accuracy of 10 cm with tha 0RAV5AT spacacraft. 
Spac if ical ly» this study will: 

(1) Davalop a data basa using available Saasat-1 altimatar data 
from 400 rapaat ravolutions. 

(2) Perform spectral analyses on all revs to assess tha 
ra idual height energy at short wavelengths. 

(3) Oanarata a global contour map of residual gaoid power using 
results from spectral analysis. 

^.aiirfliEiL 

A thorough spectral analysis of 5easat-I altimeter data on a 
global basis for suitable arc lengths serves to better delineate 
the geographic boundaries of regions with low residual gaoid 
power at wavelengths too short to be represented in CRAVSAT 
measure.ments. By spectral analysis of Seasat altimeter data at 
frequent intervals along the 3 day repeating orbit arcs. it 
should be possible to define a dense grid of residual geoid power 
(gaoid roughness) estimates for the world's oceans. 

Summary af Work Performed 

Seasat altimeter data provide a detailed and accurate 
measure of the instantaneous sea surface height ever^ 3 km elong 
the orbit track. These data comprise about 1900 orbit 

revolutions! most important of which are revolutions 1020 to 
1901f the 3-day repeat orbit. In this or^i*;, « repeat track of 
altimeter data is generated every 3 days within a kilometer or 
two cf previous tracks. Each colinear track measures essentially 
Che same geoid# but different temporal ocean surface phenomena. 
The grid of 3-day repeat tracks has a sp'.cing of 8.372 degroes of 
longitude. Brown# et al. C19813 have analyzed the variation of 
geuid roughners from these data in 9000 km arcs# finding a range 
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of values from 170 cm RM3 in the western Pacific to less than 7 
cm RMS in the southeastern Pacific. However* these relatively 
long arcs prevent accurate delineation of local areas and lead to 
considerable anisotropy between ascev^ding and descending passes 
over the same area. 

However* if we break up these 3-»fay data arcs into 2 minute 
time intervals <800 km segments)* then we can estimate the 
variability of the altimetric sea height power spectra at a 
resolution consistent with the orbit track spacing. Subsequent 
separation of the oceanic or temporal variation spectra from the 
geoid spectra would permit construction of a global mao of the 
spatial variation of geoid roughness at the same resolution. The 
span of Seasat orbits from 1199 through 1367 was found to yield 
the best redundant coverage of the marine geoid. Mhile the use 
of a 19 day span of data gives 9-fold redundancy in coverage of 
some areas* the presence of data gaps and otherwise unusable data 
in some of the arcs make the 19 day span a necessity. In fact* 
even with this long an arc* there are still soma gaps in the 
North Atlantic coverage. 


Defore useful spectral analyses can be performed on these 
data* they must be reduced to sev« surface heights relative to the 
reference ellipsoid. If long-wavelength orbit error is present 
in the altimetric geoid data* it can distort the shape and slope 
of the derived power spectrum* and cause erroneous estimation of 
the total short wavelength variance for that segment. Thus an 
accurate estimate of satellite altitude above the reference 
ellipsoid is required* a quantity which is usually not known to 
better than several meters. To overcome this difficulty* we have 
fit the sea surface heights derived from the altimeter data to a 
reference surface which approximates the shape of the geoid 
surface. This data processing is accomplished by fitting a low 
order function to the data* (see Brown C17823) over a 6000 km 
pass comprising 8 short segments. Subsequent subtraction of this 
low order function from the data results in sea surface heights 
which are a least squares best fit to the OEM-lOB geoid model 
(Lerch et al. * 1978). 


This orbit correction process is performed for 904 passes in 
the 19-day arc. Each of these passes has precisely defined start 
and stop points to ensure exact coincidence with passes from 
subsequent 3-day arcs. Passes' which had suspiciously large 
residuals 09 m RMS) relative to the GEM-lOB geoid were deleted. 
There were no more than 17 such passes out of 904 processed in 
the 19 day arc. These were caused primarily by erroneous data on 
the altimeter tapes. Each 6000 km pass was then divided into 
exactly 8 segments at prescribed points for subsequent spectral 
analy sis. 

Tl^e altimetric geoid height data in each 2 minute segment 
was tr^fated as a conventional time series for spectral analysis. 
Using « window which minimizes spectral leakage effects* a 
conventional periodogram analysis was used to construct the power 
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spectrum for wavelengths betwetn 800 km and 7 km. The power 
represented here for a given frequency is the sum of squares of 
the coefficients of the sine an<^ cosine terms of that frequency* 
and thus corresponds to a peak power. The average power* i. e. 
the mean square of the amplitudes over the segment* is just half 
the peak power. Note that the frequency is specified in cycles 
per revolution about the earth* which corresponds to the degree 
index of a spherical harmonic expansion, A power-law 
representat ion is then fitted to the spectrum between 90 and 1000 
cycles per revolution* corresponding to the wavelength range QOO 
km to 40 km. This frequency range includes not only the geoid 
signal of interest* but the mesoscale temporal variations of the 
ocean and part of the white noise floor of the altimeter 
measurement as well. Note that the wavelengths of ocean tides 
are generally longer than 800 km* so they may be considered not 
to affect this spectral analysis. 

Since the temporal phenomena and the noise tend to dominate 
the higher frequency range* we give higher weight to the lower 
frequencies in the power law fit. 8uch a weighting also tends to 
compensate for the over-representation of the higher frequencies 
in the periodogram analysis. A weight inversely proportional to 
frequency squared was adopted in fitting the power law. TtiO 
slope of the power law agrees closely with the -3.04 slope 
predicted by Kaula's rule of thumb for a linear spectrum (Nagner* 
1979). 

The power-law for each short arc must be scaled-down by its 
proportion of a complete orbit to obtain a spectral energy 
estimate comparable with a global spectrum. Ne express this as 
the fraction of the segment length in seconds compared to the 
orbital period of 6049.2 seconds* and scale the power law 
accordingly (e. g. Nagner* 1979). The scaled power law is 
extrapolated and integrated for all wavelengths smaller than 220 
km to obtain an estimate of the variance of geoid height in this 
shan't wavelength range. The square root of this quantity is 
termc'd the geoid roughness for that particular segment. This 
same analysis is performed for each of 9270 individual segments 
of the 19 day arc. Ignoring those segments which had data gaps 
greater than 6 seconds cumulative* 3011 roughness values were 
gridded (averaged) on a global map at 9 degree intervals* and 
contoured. The resulting geoid roughness map is shown in Figure 
3 of the attached manuscript “Roughness of the Marina Oeoid From 
Seasat Altimetry'*. Note that these roughness values are in terms 
of the square root of peak powers* and must be scaled by a factor 
of 0.707 to obtain roughness in terms of RMS of geoid heights. 

The attached manuscript has been submitted for publication 
in the special Seasat issue of JOR. A paper of the same title 
was presented at the AOU meeting in San Francisco. It was well 
received* and seemed to be clearly understood by the audience. A 
letter was subsequently received from Professor Rapp which took 
issue with statements in our paper regarding treatment of ORAVSAT 
truncation errors. We responded to his letter* pointing out that 
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u/e had adopced a conservative approach to the estimate of 
aliasing due to truncatloif uihcreas he had optimistically 
neglected it entirely. Professor Rapp ac knouledged that the 
aliasing error Sue to truncation had been neglectedi but uiould be 
treated in a forthcoming publication by Colombo. Our analysis in 
the geoid roughness paper had made the very conservative 
assumption that this error would be so large as to prevent 
recovery of any gravity features whose wavelength in smaller than 
about 200 km. 


In preparation for future analyses (see Recommendations > i 
the tapes of the SYNBAP5 bathymetry data base were copied at 
05FC. Associated documentation is being studied to transform 
meridional parts and Marsden sq,uares to latitude and longitude 
for correlation analyses. 

S9n^lva.La.!vg. 


220 


The variability of geoid power 
km (geoid roughness)/ has been 


between latitudes 72°N and 72°S. A 


for wavelengths shorter than 
mapped for the world's oceans 
spectral analysis of Seasat 
altimeter data# reduced to sea surface heights* has been 
performed at 2 minute intervals for 19 consecutive days of the 
3~day repeat orbit. The geoid roughness represented by these 
spectra for wavelengths shorter than about 220 km is separated 
from the total sea height variance and is displayed in the form 
of a global contour map. The global average guoid roughness is 
32 cm RliSi varying from a high in excess of 2 m RMS near deep 
ocean trenches to a low of 2 cm RMS in the southeast Pacific near 
the East Pacific Rise. This average value agrees well with 
previous estimates based on gravimetry and GEOS-3 altimetry. In 
general/ the smoothest areas in the marine geoid overlie 
relatively young seafloor adjacent mid-ocean spreading centers# 
where.' even short-wavelength topographic variations tend to be 
i sostat leal ly compensated. 


The geoid roughness map 
the error of omission in a 
In certain rather large area 
of omission should be sma 
aseismic ridges. Insofar as 
trusted# we feel that the 
GRAVSAT mission can be satis 


can also be interpreted in terms of 
satellite gravity (GRAVSAT) mission, 
s of the southeast Pacific# the error 
Her than that over trenches and 
the absolute roughness values can be 
10 cm geoid accuracy goal of the 
fled best in the southeast Pacific. 


There are additional analyses which can be performed on 
these data which are even more interesting for understanding the 
inner workings of the earth. Given the availability of detailed 
bathymetry from uhich depths along the satellite track c«m be 
inferred# it is feasible to calculate the cross-spectrum between 
the altimetric sea height and the bathymetry# and evaluate the 
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depth of compensation at various points 
This could give rise to a related 
variability of the depth of isostatic 
Implications for the global patterns 
heat flou. 


along the data track, 
map shoeing the spatial 
compensation^ ujhich has 
of crustal thickness and 


It is also possible to learn something more from the data 
processing already completed. I refer to the parameters of the 
pobiar*'lauf fit to the sea surface height covariance spectra. In 
particular/ the exponent of the fit/ ufhich varies from -2.07 near 
the East-Pacific rise to <-9.0 in the Indian ocean geoid louji can 
be related to the depth of the underlying density structure. 
This is essentially due to the more rapid attenuation of short 
uiavelength geopotential anomalies uith depth. For example/ a 
white noise distribution of density anomalies at the earth's 
surface that would exhibit a flat spectrum of geoid heights (zero 
slope)> shows a spectrum of increasingly steeper slope as the 
density distribution is positioned at larger depths. 


Some depth analyses have been performed using spectra of 
relatively long wavelength spherical harmonic geopotential 
models. Cook C19633 has shown from analysis of degree variance 
spectra of the low order geopotential that only anomalies of 
wavelength greater than 4000 km can originate as deep as the 
core-mantle interface. A higher degree model (8A0 II) was 
analyzed by Allan C17723/ who concluded that anomalies of 
wavelength greater than 7000 km are due to density fluctuations 
in the lower mantle/ while anomalies of wavelength shorter than 
7000 km can be generated by fluctuations in the asthenospher e. 
The nature of these spherical harmonic geopotential models is 
such that it is only possible to deduce density structure at 
asthenospher Ic depths and deeper/ and then only in a global 
average sense. 

In the field of aeromagnetic data interpretation/ however/ a 
different situation prevails. The existence of highly detailed 
aeromagnetic profiles and maps permits spectral analyses at 
shorter wavelengths/ and the resolution of shallower depths for 
definite areas. At the same time# the existence of distinct and 
shallow horizons such as the Curie point and the boundary between 
crystalline and sedimentary rocks/ result in relatively shallow 
depths of origin for magnetic anomalies. By postulating a 
statistical model of magnetic blocks of varying depth/ width/ 
thickness/ and magnet i zat ion/ it can be shown that the logarithm 
of the magnetic field energy spectrum is proportional to the 
depth of the ensemble of blocks. This technique has been widely 
used for determination of depth to magnetic basement. In 
principal/ this method is equally applicable to analysis of geoid 
height profiles. 

While density anomalies are not limited to such shallow 
depths/ and they do not have such distinct horizons as magnetic 
anomalies/ it should be possible to learn something about the 
variation of density structure within the crust simply by 
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analysis of the already existing set of spectral slope values. 
As uiith any neuf technique* some experimenting must be done to 
determine the proper interpretation of the derived depths. 
Houjever/ at t'le very least a global map of relative depth to 
density structure can be constructed uhich can have great value 
for interpretation of the tectonic provinces of the earth. In 
brief* uie seem to have discovered a very fertile area for 
geophysical research: one uhich promises to bear fruit very 

quickly in terms of global marine maps of compensation depths and 
depths of density structure. 
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ROUGHNESS OF THE MARINE OEOID FROM SEASAT ALTIMETRY 
R. D. Brouini W. D. Kahn*# D. C. McAdoo*# and W. E. Himwich 

Phoenix Corporation 
1700 Old Meadow Road 
McLean# Virginia 22102 

The geographical variability of short wavelength geoid power 
spectra (geoid roughness}# has been mapped for the world's oceans 
between latitudes 72°N and 72° S. A spectral analysis of Seasat 
altimeter data# reduced to sea surface heights# has been 
performed at 2 minute intervals for 15 consecutive days of the 
3~day repeat orbit. The geoid roughness represented by these 
spectra for wavelengths shorter than about 220 km is separated 
from the total sea height variance and is displayed in the form 
of a global contour map. The global average geoid roughness is 
32 cm RM8# varying from a high in excess of 2 m RMS near deep 
ocean trenches to a low of 2 cm RMS in the southeast Pacific near 
the East Pacific Rise. This average value agrees well with 
previous estimates based on gravimetry and GEOS-3 altimetry. In 
general# the smoothest areas in the marine geoid overlie 
relatively young seafloor adjacent mid*-ocean spreading centers# 
where even short-wavelength topographic variations tend to be 
isostat ical ly compensated. 


This study was sponsored by NASA Goddard Space Flight Center 
Oeodynamics Branch under contract NAS5-26S79. 


* NASA Goddard Space Flight Center# Greenbelt# Md. 
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The advent of satellite altimetry has provided an 
opportunity to map the n .rine geoid at scales as small as a few 
tens of kilometers with sub-decimeter precision. Since the fine 
structure of the geoid surface is sensitive to the depth of 
compensation of bathymetric featuresi me can begin to infer 
variations in depth of compensation and general structure of 
oceanic crust from satellite altimetry. Spectral anal^isis 
provides a convenient tool for evaluating these reiat ionsh ip s/ 
since it enables some separation of the geoid signal spectra from 
that of the altimeter noise and oceanograph ic effects. Wagner 
and Colombo C1979!I have shomn that it is possible to obtain 
accurate local estimates of the geopotential pouter spectrum from 
spectral analysis of 0EO3-3 altimeter data. Here ue apply 
spectral analysis to short ^rcs of Seasat altimeter to derive a 
preliminary map of the variability of the cumulative geoid pouter 
for scales less than 220 km 'geoid roughness). Seasat altimeter 
data are much more precise than GEOS-3 data> permitting a more 
accurate separation of the short utavelength geoid signal from the 
altimeter noise and oceanograph ic effects. The follouting 
sections detail the procedure by uthich Seasat altimeter data are 
processed to construct the geoid roughness map. This map> 
derived from 15 consecutive days of Seasat data< is interpreted 
for error of omission of satellite gravity (QRAV/SAT) missions as 
uiell as for depth of isostatic compensation of th«!; ocean crust. 
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Seasat altimeter data provide a detailed and accurate 
measure of the instantaneous sea surface height everg 3 km along 
the orbit track. These dat^ji comprise about 1500 orbit 

revolutions^ most important of uihich are revolutions 1020 to 

1S01> the 3-day repeat orbit (see Figure 1). In this orbits a 

repeat track of altimeter data is generated every 3 days within a 

kilometer or two of previous tracks. Each colinear track 

measures essentially the same geoid, but different temporal ocean 
surface phenomena. The grid of 3-day repeat tracks has a spacing 
of 0.372 degrees of longitude. Broujn/ et al, C19813 have 

analyzed the variation of geoid roughness from these data in 5000 
km arcs« finding a range of values from 170 cm RMS in the western 
Pacific to less than 7 cm RMS in the southeastern Pacific. 
However^ these relatively long arcs prevent accurate delineation 
of local areas and lead to considerable anisotropy between 
ascending and descending passes over the same area. 

Houeverf if we break up these 3-day data arcs into 2 minute 
time intervals (800 km segments)# then we can estimate the 
variability of the altimetric sea height power spectra at a 
resolution consisirent with the orbit track spacing. Subseq^uent 
separation of the oceanic or temporal variation spectra from the 
geoid spectra would permit construction of a global map of the 
spatial variation of geoid roughness at tbe same resolution. The 
span of Seasat orbits from 1155 through 1369 was found to yield 
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the beet redundant coverage of the marine geoid. The coverage 
from this span is shouin in Figure li luhich also indicates the 
size of the 800 km segments. Wh^le the use of a 15 day span of 
data gives “i“fold redundancy in coverage of some areas# the 
presence of data gaps and otheruiise unusable data in some of the 
arcs make the 15 day span a necessity. In fact# even with this 
long an arc# there are still some gaps in the North Atlantic 
coverage. 


Orbit Enror Correction 

Befare useful spectral analyses can be performed on these 
data# they must be reduced to sea surface heights relative to the 
reference ellipsoid. If long-wavelength orbit error is present 
in the altimetric geoid data# it can distort the shape and slope 
of the derived power spectrum# and cause erroneous estimation of 
the total short wavelength variance for that segment. Thus an 
accurate estimate of satellite altitude above the reference 
ellipsoid is required# a quantity which is usually not known to 
better than several meters. To overcame this difficulty# we h^.ve 
fit the sea surface heights derived from the sltimeter data to a 
reference surface which approximates the shape of the geoid 
surface. This data processing is accomplished by fitting a low 
order function to the data# (see Brown C19823) over a 6000 km 
pass comprising S short segments. Figure 1 shows an example of 
two such passes in the Pacific# their start and stop points 
denoted by the arrowheads on the orbit subtrack. Subsequent 
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subtraction of this louj ordtr function from tho data results in 
sea surface heights uhich are a least squares best fit to the 
OEM-lOB geoid model <Lerch et al.il970). 

This orbit correction process is performed for 904 passes in 
the IS-day arc. Each of these passes has precisely defined start 
and stop points to ensure exact coincidence with passes from 
subsequent 3‘^day arcs. Passes which had suspiciously large 
residuals 05 m RMS) relative to the GEM-lOB geoid were deleted. 
There were no more than 17 luch passes out of 904 processed in 
the 15 day arc. These were caused primarily by erroneous data on 
the altimeter tapes. Each 6000 km pass was then divided into 
exactly Q segments at prescribed points for subsequent spectral 
analysis. 

Spectral Ana lusts 

The altimetric geoid height data in each 2 minute segment 
was treated as a conventional time series for spectral analysis. 
Using a window which minimizes spectral leakage effects/ a 
conventional periodogram analysis was used to construct the power 
srectrum for wavelengths between 800 km and 7 km. The spectrum 
for a typical segment is shown in Figure 2. The power 
represented here for a given frequency is the sum of squares of 
the coefficients of the sine and cosine terms of that frequency/ 
and thus corresponds to a peak power. The average power/ i. e. 
the mean square of the amplitudes over the segment/ is just half 
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the peak power. Note that the freq.uency is specified in cycles 
p«r revolution about the earth/ which corresponds to the degree 
index of a spherical harmonic expansion, A power-law 
representation is then fitted to the spectrum between 50 and 1000 
cycles per revolution# corresponding to the wavelength range 800 
km to 40 km. This freq.uency range includes not only the geoid 
signal of interest# but the mesoscale temporal variations of the 
ocean and part of the white noise floor of the altimeter 
measurement as well. Note that the wavelengths of ocean tides 
are generally longer than 800 km# so they may be considered not 
to affect this spectral analysis. 


Since the temporal phenomena and the noise tend to dominate 
the higher frequency range# we give higher weight to the lower 
frequencies in the power law fit. Such a weighting also tends to 
compensate for the over-representation of the higher frequencies 
in the periodogram analysis. A weight inversely proportional to 
frequency squared was adopted in fitting the power law. The 
resulting power law# and its parameters for a typical segment in 
the South Atlantic# are shown in Figure 2. The slope of the 
power law agrees closely with the -3.04 slope predicted by 
Kaula's rule of thumb for a linear spectrum (Wagner# 1979). 
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The po(uer-*laui for each short arc must be scaled~douin by its 
proportion of a complete orbit to obtain a spectral energy 
estimate comparable uiith a global spectrum. Ue express this as 
the fraction of the segment length in seconds compared to the 


orbital period 

of 6045. 2 seconds* 

and 

scale the power 

law 

accordingly (e. g 

Wagner* 1979). 

The 

scaled power law 

is 

extrapolated and 

integrated for all 

wavelengths smaller than 

220 

km to obtain an 

estimate of the variance 

of geoid height in 

this 

short wavelength 

range. The square 

root 

of this quantity 

is 


termed the geoid roughness for that particular segment. This 
same «ru.£!^iis is performed for each of 5270 individual segments 
of the? IS day arc. Ignoring those segments uihich had data gaps 
greater than 6 seconds cumulativei 3011 roughness values mere 
gridded (averaged) on a global map at 5 degree intervalsi and 
contoured. The resulting geoid roughness map is shown in Figure 
3. Note that these roughness values are in terms of the square 
root of peak po(uers< and must be scaled by a factor of 0. 707 to 
obtain roughness in terms of RMS of geoid heights. 


There are several assumptions 
the geoid roughness which deserve 
the proper description of the 
Uhile this has not been verified* 
it is of the same form as Kaula's 
and is easily integrable. Until 


involved in the derivation of 
examination. The first is that 
geoid spectrum js a power law. 
a power law was adopted because 
rule for the l(fng wavelengths 
we can perform a more complete 
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tt»par«tion of th® 8p»i:tra of th« octan surfac® variatione from 
tha gooid spectra^ tha adoption of a simplo pouar lain seoma an 
adaq.uata choica. 

The second aasumption la that the inverse fraq.uencg sq.uared 
uiaighting fairly attenuates the masoscale ocean surface variation 
spectra (SO to SOO km yavelangth) and the altimeter noise spectra 
uhilo allouiing an accurate representation of the short utavelength 
geoid spectra. This is probably over-c imp 1 ist ici but it is easy 
to implemantf and it gives pouter laut parameters uthich agree utelli 
on average* uith longer utavelength spectra. Given the precision 
of the pouter spectrum* it is an appropriate rule for a 
preliminary mapping of geoid roughness. 

Finally* uie must consider the uncertainty of the spectral 
analysis technique itself. It is utelX knouin that the use of any 
spectral analysis technique requires certain assumptions* o. g. 
periodicity of the data in the time series. Each technique also 
has its oum set of heuristic rules which must be applied to 
obtain reason«ible results. For example* the Blackman and Tukey 
C19983 method requires pre-whitening and smoothing windows. 
Leas t-squares spectral analysis* per iodograms* fast Fourier 
transforms* etc» have similar special assumptions. If all these 
techniques were applied to the same data segment* the variation 
in results would probably exceed any error due to the above 
assumptions. Th« most important factor in conducting an 
investigation uving spectral analyses is consistancy and 
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simplicity. Insofar as data are processed uniformly and 
minimally/ results should be valid relatively/ if not absolutely. 
Therefore/ uhile the absolute values of the gcoid roughness 
generated by this process may be in error/ their relation to 
their neighbors (geographic variation) remains valid. 

Risykia. 

The geoid roughness map of Figure 3 is the spatial average 
of the roughness values of ascending and descending segments/ 
redundant in general by a factor of 5. This map is displayed on 
a Van der Orinten projection for comparison with the physiography 
of National Geographic Society maps. The geoid roughness values 
are color coded (in cm) according to the key in the upper right 
corner of Figure 3. 

The largest contrast in this map occurs between the geoid 
roughness of smooth regions like the southeast Pacific and rough 
areas like trenches. It is interesting to examine in greater 
detail the spectra of arc segments in these two regions. Figure 
4 shows the spectra of a typical trench segment/ superposed on 
spectra of a segment over the smooth flank of the East Pacific 
rise. The essential difference in these spectra is their slope/ 
with trench regions exhibiting larger slopes (in an algebraic 
sense) than the global average. Figures 5a and 5b show the 
residual altimetric sea heights for the selected trench segment 
and the segment near the East Pacific rise. 
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The pouier lauj fit and integration shown in Figure 2 has been 
performed for each of 3011 arc segments. The global average 
slope of these power laws is ~3. 409« and the global average 
roughness is 40 cm (32 cm Rt1S>« yeilding a global average power 
law given by 


-3. 409 2 -1 

P *= 132195 f (m Ecycles per rev3 > 


( 1 ) 


This power law« properly scaled for comparison with short arc 
spectra> is also plotted in Figure 4. The integrated power for 
wavelengths shorter than 220 km has been computed for each of 


wif«T9S. I II * a vw A u t«i i iiaa wcrv;fi 


ll*»Ayil«r\« VW WIII7 


mid-point of the short arc segment and gridded in a 5 degree by 5 
degree lat i tude-longitude grid. Thus 5 degree average variances 
are formed from the geoid variances of both ascending and 
descending arc segments/ with up to five-fold redundancy in some 
areas. The variation of roughness between colinear segments is 
about 3 cm. Except near long linear features such as the 
Tonga-Kermadec trench/ the variation between ascending and 
descending segments is remarkably small. The most severe 
anisotropy/ 150 cm vs. 40 cm/ occurs for this trench which 
parallels the descending arc segment. In other less lineated 


areas/ the anisotropy is less than 10 cm. 
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Thff short wavelength geoid roughness estimated in this study 
must of necessity be due to near-surface density variations of 
short wavelength in the crust and lithosphere. For example* we 
cite the spatial correlation between the distribution of geoid 
roughness mapped in Figure 3 and the degree of roughness of 
seafloor topography inferred from a global bathymetric chart. 
Deep ocean trenches represent sites of substantial short 
wavelength bathymetric relief and appear as highs on the geoid 
roughness map* Figure 3. This correlation is expected in view of 
the correlation generally found (McKenzie and Bowin C19763) 
between bathymetry and shipborne gravimetric observations* and 
between bathymetry and altimetric geoid observations (e. g. 
Leitao and McOoogan C19743* Wagner C19793* Chapman and Taiwan! 
£19793). 


This geoid roughness map* however* contains more than just a 
confirmation of the gravi ty/bathymetry correlation. Regional 
variations in the depth of isostatic compensation can also be 


inferred 
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from the map. Over regions uhere isostatic compensation occurs 
only at substantial depth/ if at all* high bathymetric relief 
eq.uates to large values of short wavelength geoid variance. 
Conversely over regions where compensation depths are relatively 
shallowi high bathymetric relief does riot result in large values 
of geoid roughness. For example/ the East Pacific rise crest/ 
while relatively rough in terms of bathymetric relief/ is one of 
the smoothest areas of the geoid/ due in part to its shallow 
isostatic compensation. This qualitative observation is 
quantified in the following paragraphs. 



In order to estimate the effect of compensation depth on the 
geoid power spectrum/ let us assume that bathymetry varies in one 
dimension only and has a power spectrum equivalent to that of 
band-limited white noise (cf. Blackman and Tukey C19583 p79). 
The power spectrum is given fay 




where f 

is frequency 
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effective cut-off frequency which is much lower than 
bathymetric power spectrum may be approximated by 
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(3) 


In order to obtain the output gsoid signal/ the input bathymetric 
signal must be passed through a linear system (or filter) which 
is represented by a frequency response (or transfer) function/ 
Y(f). Then the p ower spectrum/ the geoid signal is 
given by 


Pg(f) « lY(f)|2p^(f) 


(4) 


The transfer function/ Y(f)/ in the case of uncompensated 
topography may be written as 

Y^(f) « Cexp(~2Tr|£ |d)[|f I ]"^ O) 

where d is ocean depth and C is a constant described below. This 
function is obtained by multiplying the function defined in 
equation (12) of McKenzie and Bowin C19761/ by the transfer 
function 

S(f) = [2Trglf|]"^ 


between gravity and the geoid (see Chapman C19773). To represent 
the case of isostatic compensation/ the transfer function 
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(7) 


is used/ where T is the conpensation depth below seafloor and 


C » ApG/g 


( 8 ) 


where Ap is crustal density minus water density. The nodel in 
eq.uation 7 is equivalent to that of simplified Airy compensation 
given by equation 13 of McKenzie and Bowin C19763 or equation 3 
of Chapman and Taiwan! C1979D. The total geoid energy (or 
roughness) residing at frequencies above lf| 


f is 
o 


, A 


(Odf 


(9) 


where Pg(f) is given by equation 4 after substitutions from 
equations 3 and 5 or 7 as appropriate. The form of the integral 
in equation 9 is the exponential integral £ 2 ^^^ defined and 
evaluated in Abramowitz and Stegun C1970]/ where z=4Ttf^^d in the 
uncompensated case. 


The square root of the short wavelength geoid roughness for 
various depths of compensation/ T/ is shown in Figure 6 . This 
power function has been normalized by the residual geoid 
roughness from the uncompensated case. Values of d=5 km and 
fo=0.0045 cy(.les/km have been adopted. The square root is 
displayed to allow a direct comparison with the geoid roughness 


map of Figure 3. This model suggests that short wavelength geoid 
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roughness mill be subs 
are 30 km or less, 
deeper than 40 km» the 
of the uncompensated 
spectra are not uihite< 
(uhich are more lilceiy 
the reduction of geoid 
(smaller than reality) 


tantially reduced if compensation depths 
□n the other hand< if compensation occurs 
geoid roughness uill essentially be that 
case. Insofar as actual bathymetric pouer 
but are stronger at long luavelengths 
to be compensated > f the model estimates of 
roughness by compensation are conservative 


D.g.efl....P.£aan .J.r. ?nc .h.f s 


Trenches are generally thought to be relatively 
uncompensated features ( Vening-Meinesz C19641). Certainly their 
effective depth of compensation is relatively large* if not 
infinite. From Figure 6 ue see that the geoid roughness in this 
case is not appreciably reduced by isostasy. Furthermore* 
trenches generally constitute exceptionally high bathymetric 
relief. In this case it is not surprising that the major maxima 
in the geoid roughness map occur over large ocean trenches. In 
fact* the global maximum of geoid roughness (more than 200 cm ) 
occurs over the Challenger Deep segment of the southern Marianas 
trench* uhich is the deepest point in the oceans. In general* a 
semi-quanti tativ« relation is observed between the maximum depth 
of a given trench «nd its geoid roughness. 


Page 
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Another feature wf the geoid roughness map is the broad 
region of loui roughness <<20 cm ) associated uith mid-ocean 
ridges. Consideri for exampl^e* the East Pacific Rise and its 
ridge flanks. Cochran C19793 analyzed gravity and bathymetry in 
this region and found that isostatic compensation occurs at quite 
shalloui depths* even shallower than those he estimated for the 
Mid-Atlantic ridge. In the case of the simplified Airy model* 
Cochran shows that a compensation depth of 10-20 km best fits the 
data for the East Pacific. According to model results in Figure 
6t isostatic compensation at depths of 10-20 km should 
effectively halve the residual geoid roughness* compared to that 
of similar uncompensated bathymetric relief. 

The region encompassing the Mid-Atlantic ridge and its 
flanks shows moderate values of geoid roughness <10-50 cm in the 
north and 10-30 cm in the south). The seafloor topography does 
not seem to be substantially rougher than thjit of the East 
Pacific Rise* but Cochran C19793 finds simplified Airy 
compensation depths of 20 to 50 km. Thus it appears that the 
disparity between geoid roughness over the Mid-Atlantic ridge and 
the East Pacific Rise is due in part to different depths of 
compensation. Cochrane’s analysis does not distiguish between the 
northern and the southern Mid-Atlantic ridge* but in view of the 
smoother geoid in the south it is interestinj} to note that the 
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southern ridge is generally believed to spread at a faster rate 
than the northern ridge. 



The aseismic Walvis ridge is manifest in Figure 3 as a band 
of moderate geoid roughness (30-50 cm >. Directly to the south* 
the roughness rises to 50-60 cm over the Discovery Tablemount* a 
feature quite unrelated to the Walvis ridge. Although different 
models of isostasy have been proposed for Walvis ridge (Detrick 
and Watts C1979J* Angevine and Turcotte C19S0!1>« there is generc»l 
agreement that the depth of compensation is shalloui (15-30 km^ 
for the Airy model. Figure 7 shous a sea surface height profile 
for a pair of tuo-minute pass segments crossing the Walvis ridge 
and Discovery tablemount. Also displayed is a bathymetric 
profile adopted frum Figure 4 of Detrick and Watts C19793. A 
moderate value of geoid roughness (28 cm > is obtained for the 
Walvis ridge segment* compared to 82 cm for the Discovery 
tablemount segment* even though the bathymetric relief (3 to 4 
km> is comparable. The shallow compensation depth is apparently 
responsible for the modest value of geoid roughness over Walvis 
ridge. 


In contrast* seamounts are not compensated at such shallow 
depths. Watts 119781 indicates simplified Airy compensation 
depths between 30 and 90 km for the Hawa i ian-Emporer seamount 
chain* a feature which shows geoid roughness in excess of 100 cm. 
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Bermuda also appears as a significant high in the geoid roughness 
map 080 cm )i u/hich tends to confirm the finding of Haxby and 
Turcotte C19783 of deep compensation for Bermuda. A striking 
high in the geoid roughness map OllO cm ) over the Tuamotu 
archipelago suggests that compensation is not shallou in this 
region either. 


E rror of-Omission for Satellite Qravitu Missions 


The geoid roughness map can also be interpreted in terms of 
the error of omission of a satellite gravity (GRAVSAT) mission. 
The GRAV-SAT mission objective is a detailed and accurats map of 
the earth's gravity field. This map is to have a resolution of 1 
degree (220 km)» and a preciaion corresponding to 10 cm RMS in 
geoid height. The planned configuration is for tmo satsllitesi 
seperated by about 330 km< in the same circular polar orbit at 
160 km altitude. By precise measurement of the change in 
relative velocity betu>een the satellites as they respond to 
anomalous gravity# one can calculate a global map of area mean 
gravity parameters. Uhile the satellites' response to gravity 
anomalies is theoretically limited only by the system integration 
time and sensitivity# the parameterization of these measurements 
as area means of a certain size limits the mission accuracy. One 
limit is the omission error# which is t^#)e cumulative effect of 
those gravity a^^omalies of half~waveleng th shorter than tha size 
of the area mean gravity parameters# another limit is the error 
due to aliasing cP the shorter wavelength orbital measurements 
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into the longer wavelength parameters. Our geoid roughness map 
gives a measure of the error of omission for a GRAVSAT mission 
parameterized for q.uasi~point q,uantitiesi insofar as our 
roughness values (wavelengths shorter than 2 arc degrees in 000 

o 

km arc segments) represent the geoid variance about a 1 area 
mean value determined by averaging of a dense set of uniformly 
distributed ORAVSAT measurementi. . 

To quantify the expected reduction in our roughness values 
which might accrue from the use of true 1 area averages^ 
consider the error of omission predicted for point measurements 
from Kauia's rule of thumb. This rulei integrated for all 
spherical harmonics of wavelength shorter than 220 kmi yields an 
error of omission of 36 cm RMS (see Chovitz* 1973). for these 
wavelengths. If we assume the geopotential spectra behave as 
predicted by Kaula^s rule of thumb for these short wavelengths< 
then we can apply Pellinen^s averaging operator to these point 
value statistics to obtain their equivalent in area means For 
examp le: 

<N^> = Re^Pj,^ (2il+l) a| (10) 

‘nil' 

(see Tscherning and Rappi 1974> for * 0.564 degrees (the 
spherical distance radius of an area equivalent to 1 degree by 1 
degree)i and where the P are the Legendre polynomials and <j| is 
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the mean square fully-normalized geopotontial coefficient of 
degree K, and Re is the radius of the earth. 

Evaluation of equation 10/ using Kaula's rule for <J ^ and 
summing between degrees ISO and 1000/ yields an error of omission 
in area mean geold heights of 17 cm RMS instead of the 36 cm RMS 
predicted for point measurements. A similar reduction should 
accrue for the geoid roughness values computed and mapped in this 
study. On the other hand/ our global average roughness value of 
32 cm RMS is in excellent agreement with the degree variance of 
31 cm RMS obtained by Wagner and Colombo C19793 from gravimetry 
and GEOS-3 altimetry for wavelengths shorter than 220 km. 

Never the 1 essi it appears that the error of omission for 
GRAVSAT may be less than 10 cm RMS for large areas of the oceans. 
The southeast Pacific ocean geoid is potentially the most 
accurately mapped by GRAVSAT/ and the areas near trenches# and 
aseismic ridges less 


so. 
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CONCLUSIONS 

The spectral analysis of Seasat altimeter data over short 
arc segments has been shoun to yield geophysically interesting 
results. To a degree* the geoid roughness map generated from 
these analyses can be used* together uith a bathymetry chart* to 
Ioc4*te areas of shallow isostatic compensation. 


The geoid roughness map can also be interpreted in terms of 
the error of omission in a satellite gravity (GRAVSAT) mission. 


In certain rather large areas of 
of omission should be smaller 
ascismic ridges. Insofar as the 
trusted* me feel that the 10 
GRAVSAT mission can be satisfied 


the southeast Pacific* the error 
than that over trenches and 
absolute roughness values can be 
cm geoid accuracy goal of the 
best in the s.*utheast Pacific. 
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Distribution of Seasat altimeter data tracks for 15 
days from September 15 to 30/ 197Q. 


Typical altimetric sea height power spectra 
showing the power-law fit compared to Kaula^s rule. 


Geoid Roughness Map: Global distribution of the 

sq,uare root of geoid peak power for wavelengths 
shorter than 220 km. 


Altimetric sea height profiles over areas of high 
and low geoid roughness. 


Comparison of spectra from areas 
roughness showing the difference 


of high and low geoid 
in slope. 


Depth of compensation versus geoid roughness from a 
white-noise model of bathymetry assuming simplified 
Airy compensation. 

Sea height profile and bathymetry for Discovery 
Tablemount and Ualvis ridge. Note the difference in 
geoid undulation amplitude for comparable 
bathymetric relief. 
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South Sandwich Trench 
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